Changes in mobility patterns during the COVID-19 pandemic in Zambia: Implications for the effectiveness of NPIs in Sub-Saharan Africa

The COVID-19 pandemic has impacted many facets of human behavior, including human mobility partially driven by the implementation of non-pharmaceutical interventions (NPIs) such as stay at home orders, travel restrictions, and workplace and school closures. Given the importance of human mobility in the transmission of SARS-CoV-2, there have been an increase in analyses of mobility data to understand the COVID-19 pandemic to date. However, despite an abundance of these analyses, few have focused on Sub-Saharan Africa (SSA). Here, we use mobile phone calling data to provide a spatially refined analysis of sub-national human mobility patterns during the COVID-19 pandemic from March 2020-July 2021 in Zambia using transmission and mobility models. Overall, among highly trafficked intra-province routes, mobility decreased up to 52% during the time of the strictest NPIs (March-May 2020) compared to baseline. However, despite dips in mobility during the first wave of COVID-19 cases, mobility returned to baseline levels and did not drop again suggesting COVID-19 cases did not influence mobility in subsequent waves.


Introduction
The COVID-19 pandemic has resulted in an unprecedented global public health crisis that has drastically changed human behavior, particularly the travel patterns of individuals.At the start of the pandemic in early 2020, many non-pharmaceutical interventions (NPIs) were put in place to minimize the spread of SARS-CoV-2 [1].Some NPIs including travel restrictions, business/school closures, and social distancing were implemented to varying degrees globally and have been shown to minimize the likelihood of onward transmission [2].Spurred by the public health relevance and importance of predicting the likely trajectory of SARS-CoV-2, datasets quantifying human mobility have also become more readily available [3].This has led to a deluge of analyses focused on quantifying reductions in human mobility and implications for evaluating the effectiveness of NPIs [4].However, despite global efforts, to date there have been few analyses focused on Sub-Saharan Africa (SSA).The pandemic has shown substantially different patterns in SSA compared to other regions with overall lower reported cases and mortality rates despite serological results suggesting overall high rates of infection [5].Limited surveillance may be a possible cause for the differences in the pandemic to date, however other factors such as protection from endemic coronaviruses, an overall younger average age resulting in less severe infections, and the early implementation of NPIs relative to the number of reported cases have been posited as possible reasons for some of the observed differences [6].Disentangling these factors requires detailed data on human mobility and behavioral patterns from countries in SSA.
Understanding how human behavior has changed during NPIs in SSA has been explored, but primarily via aggregated data as a proxy for mobility such as Google Mobility Reports [7].Given biases in who is represented in these data based on Google app usage, these data may be presenting a biased view of behavior.Kishore et al. (2021) point out several of these biases in the limitations of their study [8].Further, these data often have spatial data gaps or are aggregated to represent mobility throughout the entire country, which may hinder detailed analyses highlighting sub-national patterns of transmission.In addition, spatial heterogeneity in transmission is further impacted by limited testing availability which limits approaches to understand transmission dynamics directly from confirmed case or death data.As of early 2022, vaccination in Africa remains limited with just under 16% of persons fully vaccinated, compared to some high-income settings where fully vaccinated proportions are greater than 60% [9][10][11].In these settings, travel and social distancing restrictions may be more readily implemented to further stem the spread of SARS-CoV-2.Additionally, changes in the virus and the emergence of new variants alter our landscape of protection and susceptibility requiring NPIs to be implemented despite widespread vaccination [12].
Using detailed data on mobility patterns quantified from mobile phone calling data, we provide the first spatially resolved analysis of mobility during the COVID-19 pandemic in Zambia from 01 March 2020 to 07 July 2021.We couple this data with information on reported cases andimplementation of various NPIs to build a model of mobility and a transmission model to understand where, how, and when travel has changed sub-nationally since March 2020.We further investigate implications for the effectiveness of travel restrictions moving forward as the course of the pandemic continues to spread through populations and prolong morbidity and mortality.

Data on mobility
Deidentified, anonymized, call detail records (CDRs) were obtained with permission from a leading mobile phone operator in Zambia.CDRs from 108 districts across all 10 Zambia provinces were included in the final data (Fig 1A).Individual subscribers were assigned a most used mobile phone tower per day based on the routing tower where most of their records were recorded.Trips were defined based on the most used mobile phone tower on subsequent days where a trip was counted if these two towers were different.If these two most used towers were the same, then this was recorded as a "stay" where a subscriber did not change their location.CDR data was processed further by aggregating individual movements both spatially by district and temporally by day.Aggregation was conducted at the district and province level to construct origin-destination matrices of movement between districts and provinces.Overall, between 566,685-871,851 unique mobile phone subscriber IDs were recorded per day from 1 March 2020 to 7 July 2021.Google Community Mobility Reports were obtained via open access from Google at https://www.google.com/covid19/mobility/[13].The output from Google's data focuses on national percent changes from baseline and is aggregated into several categories including changes in visits retail stores, grocery stores, parks, transit, workplace, as well as residential mobility.To compare with the mobile phone data, we only looked at the changes in workplace mobility as a proxy for overall national mobility.

COVID-19 data
Daily case data for Zambia were obtained from the COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (https://github.com/CSSEGISandData/COVID-19).Case data were available by district from 18 March 2020 to 22 February 2021 (Fig 1B) [14].

Non-pharmaceutical interventions (NPIs)
Data on dates and duration of NPIs put into place during 2020 were obtained from the Oxford Covid-19 Government Response Tracker [1].Initial NPIs were implemented nationally in Zambia from March 2020 to May 2020.The main restrictions included workplace and school closures, stay at home orders, and limitations on gatherings.The stringency level of NPIs in Zambia were the strictest between March and May 2020, but later varied throughout the course of the pandemic.NPIs implemented from March to December 2020 can be defined as: For some analyses, each wave was then averaged to obtain a district-level mean number of trips during that time period.Since 2019 mobile phone data was not available, the time period from 01 Sep-31 Dec 2020, where mobility was stable and no NPIs were in place, was used as a proxy for "baseline" mobility.While this time period does not reflect a 'true' pre-pandemic baseline value, given the time period of data we were able to access, we used the most stable period of mobility during this time when no restrictions were in place to help ground our analysis and compare across COVID-19 waves.Additionally, this period came before the next large wave of COVID-19 cases in early 2021 (Wave 3).We compared the percent change in overall mobility between Wave 1 and baseline using different baseline variations and did not find any qualitative differences between variations (S1 Table ).
Since the mobile phone data did not show any drops in mobility during Waves 2 or 3, we focused our analysis on mobility and NPIs implemented during Wave 1.For analyses comparing mobility during the different waves and baseline, the data were normalized to the baseline values.For some analyses, Welch's two-sample t-tests were used to compare groups.Lastly, we conducted a simple linear regression using public holidays and school holidays to see if there were changes in mobility.

Mobility model
To investigate the rate of decreasing mobility in the initial stage of the pandemic and the effect of explanatory variables on this decrease, we fit a Bayesian model of mobility to the observed trip counts.For each pair of districts, the rate of trips from district i to district j on day t, γ ijt , was modelled as a fraction of the baseline travel rate for this pair, r ij , The fraction q ijt represented the cumulative daily decrease in travel in the district, modified by the effect of explanatory variables, where c k > 0 was the decrease in travel on day k, X ij were explanatory variables and β was a vector of effects to be learned.Euclidean distance was used as the primary explanatory variable in this analysis.A Poisson likelihood was used for the number of trips observed each day from district i to district j, where p i was the population in district i and the model was completed by applying priors to the model parameters.Rates of decrease were parameterized on the log scale with an Autoregressive model, AR(1) prior, ϕ = 1, c = 0, and ε = 0.5.A standard Normal prior was placed on β.The model was implemented in the R programming language using the TMB package and maximum a posteriori estimates were obtained using the Laplace approximation [15].

Transmission model
To compare the impacts of different mobility information on simulations of SARS-CoV-2 transmission dynamics, we used a spatial semi-mechanistic time-series transmission model that integrated each mobility data set [16].Since short time scales were the focus of this analysis, we deemed a simple SIR model to be sufficient to draw conclusions.Key parameters in the model were M, α, and β 0 ; where M represents the mobility matrix that approximates connectivity between districts, α represents a scaling parameter to account for the discrete time steps, and the β 0 represents the transmission coefficient, which can be approximated by R 0 , which we are assuming to be 2 for this analysis (sensitivity analysis exploring different values of can β 0 be found in (S2 Fig) [17].For our model we set α to 0.99 suggesting near homogenous mixing patterns as has been done for other similar models [16,18].District-level population data were sourced from the Zambia Statistics Agency and are publicly available at https://www.zamstats.gov.zm/.Initial conditions were set to be 2% of persons infected in Lusaka district only and 95% of the persons in each district were susceptible.We investigated the impact on SARS-CoV-2 cases based on two mobility matrices informed by either the Google Mobility Reports or the mobile phone data.The transmission model in this analysis required values in the mobility matrix to fall between 0 and 1.To generate these mobility matrices, we first calculated the national average reduction of mobility in Zambia during Wave 1 compared to baseline for both Google Mobility Reports and the reduction in mobility recorded for district-pairs in the mobile phone data.Second, we multiplied these weights to the baseline mobility matrix that contains unique mobility estimates for each district pair to create a new matrix for each dataset with the Google Mobility Reports scaling travel the same nationally whereas the mobile phone data scaled travel based on the recorded reduction by district pairs.If there was no value (no trips reported), we set a value of 2 to all for no missing values.Further, since intra-district travel accounts for the majority of trips and it was unclear how to calculate this value for Google Mobility Reports, we fixed this value at 80% (sensitivity analyses in (S3 and S4 Figs).We then rescaled the inter-district trips (20% in the case where the diagonal of the origin-destination matrix was 80%).Lastly, the data were normalized with minimum-maximum scaling methods (calculated by (x-min(x)/(max(x)-min(x)), where x is the rescaled number of trips).The min-max normalization means that the re-scaling by either the Google Mobility Reports or mobile phone data does not result in the same matrix.While the same baseline mobility matrix was used for both datasets, the magnitude of reduction in inter-district mobility that was applied to the baseline matrix differed.For the Google Mobility Reports data, we used the nationallevel scaling for the overall reduction and rescaled the entire matrix based on these.For the mobile phone data, we used the actual reductions observed for each province these data.Further, we explored using population weights on data, but did not find it qualitatively changed the results and did not include any additional weighting on the data in the final model (S5 Fig).
The number of persons susceptible (S i,t ) or infected (I i,t ) at each time point, t (one time point = one week), for a district, i, can be estimated by: Where, Infected individuals recovered by the next time step in our discrete model.Each model was run fifty times and averaged for further analysis.We also estimated the time until the first 10, 50, and 100 cases were reached in each district and compared output between mobility datasets and stratified by traffic density.Number of cases were for overall infections.
Regression Analysis.To explore seasonality in movements, we fit a simple linear regression with total number of national trips per day as the outcome and national holidays and school breaks were the predictors.Since school breaks include longer periods of time between terms, they were treated as their own variable.We considered the impact of holidays on changes in mobility using a vector autoregressive model to capture temporal structure in the data.Given the strong autocorrelation (see S6 Fig) , we considered a range of lags representing 1, 7, and 12 days.A single model for a given lag was run considering both national and school holidays as binary temporal variables.National and school holidays were not significant in any of the three models (see S2 Table ).

Total number of trips per day per
Where, k indicates the number of lags (k = 1, 7, 12 for this analysis).

Ethics statement
Study approval was granted by the Johns Hopkins Bloomberg School of Public Health (IRB 15892) which deemed the work exempt.Travel within provinces also saw reductions, but to a lesser degree than inter-province (Fig 2A and 2B).Overall, intra-province travel decreased during the first wave by 23-52% of baseline (Fig 2B).The provinces where the largest number of trips were recorded saw fewer reductions in travel after restrictions were implemented.For example, travel within Lusaka, Copperbelt, and Southern Provinces saw reductions of 36%, 32%, and 23% during Wave 1, respectively.In contrast, intra-province travel within provinces with the fewest number of trips recorded saw larger reductions than those in more traveled routes (North-Western and Muchinga, experienced reductions of 42% and 47%, respectively).Intra-district mobility did not see the same reductions as inter-district mobility (reductions of 3-80% (IQR: 24-46%) during Wave 1 compared to baseline.Additionally, there were only two intra-district routes saw increases during the Wave 1 compared to baseline.This is further illustrated by the ratio of inter-to intra-province travel which fell on average 64% during the first wave (see (Fig 2C), but with differences across routes.Further, our model assumed 80% intra-district travel and, when testing the sensitivity of this assumption, we found that an increase in intra-district travel did not produce different results when looking at time to first 10 SARS-CoV-2 cases and total number of SARS-CoV-2 cases (S5 Fig) .During Wave 1, more frequently traveled routes had greater reductions in inter-province travel and lower reductions in intra-province travel.This is likely due to NPIs such as stay-at-home orders, restrictions on public events, and limitations on gatherings put into place at the beginning of the first wave, discouraging inter-province travel.

Mobility returned to pre-pandemic levels within a few months
Following the lift of the strictest NPIs on 08 May 2020, the number of trips and routes of travel returned to baseline by July 2020.Intra-province travel recovered faster than inter-province travel (intra: average = 2.5 months, IQR: 1.5-3.2months; inter: average = 3.3, IQR 1.6-4.9months) (Fig 2D).The top 10 most traveled routes recovered faster (average 2 months) whereas the 10 less frequented routes recovered slower (average 2.43 months); however, difference in recovery between the most and least commonly traveled routes was not significant (Welch's t-test p-value = 0.29).Across all routes, 15% of district routes took 6 months or greater to return to baseline.However, these routes were rarely traveled during the baseline period with an average of 1.97 trips per day and were likely not heavily trafficked during the pandemic.Once mobility returned to baseline following the lifting of NPIs during the first wave, mobility remained stable during subsequent waves, including during time periods of large national holidays with corresponding breaks in school terms (Fig 3A -3C), and large reductions in inter-and intra-province mobility were not observed.We did not find any significant relationship between school breaks or national holidays on mobility (estimates = 4.178 and -7.298; p-values = 0.0931 and 0.1343, respectively).

Mobility did not decrease in response to COVID-19 case reports
Unlike other settings globally, COVID-19 cases were less common in Zambia with few districts (25/116) reporting at least one confirmed case by the end of May 2020.Among these districts, including some with the first and highest number of cases reported (Lusaka, Lusaka Province and Nakonde, Muchinga Province) we did not see any differences in the reduction in mobility or return to baseline travel compared nationally (Welch's t-test p-value = 0.86).Further, during the subsequent waves starting in July 2020 and Jan 2021, we did not see any additional reductions in mobility.
Given variability in testing and reporting across Zambia at the start of the pandemic, we focused our subsequent analyses in Lusaka where testing was more readily available.During the period following the first case report in Lusaka on 18 March 2020 and the end of Wave 1, traffic from other provinces into Lusaka decreased, on average, 19.6% compared to baseline (Fig 4 ); however, reductions into Lusaka did not significantly differ from travel in other districts and provinces (Welch's t-test p-value = 0.08).Travel to Lusaka District from other districts within the province decreased by 67% on average compared to baseline; directly following the first case report, mobility dropped across all districts within Lusaka Province and remained low during the first wave until they eventually rebounded around May and June 2020 where it remained stable for the duration of the year (Fig 4).Further, intra-province travel in Lusaka Province had significantly greater reductions in district-level mobility compared to other provinces (p = 0.001).

Mobile phone data better illustrates district-level mobility compared to Google Mobility Reports
Google Mobility Reports have served as a popular source of data for analyses aiming to incorporate human mobility data given their global availability.However, how these data compared to mobile phone data has not been widely explored in SSA.Overall, Google Mobility Reports estimate lower magnitudes of change in mobility compared to the mobile phone data used in this analysis, however, a similar trend is consistent across datasets during Wave 1 (Fig 5A).It is important to note that Google Mobility Reports were available nationally, and thus cannot be used to understand district or province level heterogeneity (low versus high trafficked routes, for example).Additionally, the mobile phone data report larger decreases in less frequently traveled routes, which is likely driving national reductions to higher magnitudes.Google data suggests a -6.91% mean difference in mobility during Wave 1 compared to baseline whereas the mobile phone data suggests a -59.59% mean difference.Since Google Mobility Reports are national, the impacts of less traveled, subnational routes may not be contributing to overall reductions as dramatically as seen with the mobile phone data.Lastly, most Google application users are likely in more urban areas, which, despite having great reductions during Wave 1, still maintained a higher level of mobility; this prevents the data from accounting for less frequently traveled trips to rural areas.
To further illustrate the impact of using both types of data, we simulated transmission dynamics of SARS-CoV-2 during 2020 using mobility data from Google or mobile phone data.We performed a simple simulation reflective of the overall patterns of the initial stages of the pandemic in Zambia where initially only cases were reported in Lusaka district in March 2020.Further, we investigated the time until at least 10 cases were reported in all other districts either informed by Google or mobile phone data.Despite differences in the reductions in travel between the data sets, there were no large differences in time to the first 10 cases using Google data or mobile phone data (   projected earlier case counts in 47% of districts, mobile phone data projected earlier case counts in 45%, and 8% of districts had no difference between datasets used.Despite a small difference, these results align with our expectations that Google would produce earlier estimates of cases rising given that the Google data suggested smaller reductions in mobility, which would have contributed to an increased risk for introduction events in other districts.We further stratified by travel frequency to each district defined as low, moderate, or highly trafficked districts (Low: �2,240 total baseline trips, Moderate: 2,241-7,741 total baseline trips, and High: �7,740 total baseline trips).We see that as we increase the threshold of cases (time to 10 cases, time to 50 cases, and time to 100 cases), the difference in time to cases becomes less prevalent between datasets (Fig 5B -5D).

Discussion
Human mobility has been widely impacted from global to sub-national scales throughout the course of the pandemic.Dissecting the patterns of human mobility during the pandemic can help quantify the impact of NPIs on mobility as well as how changes in travel can influence SARS-CoV-2 transmission.Here, we aimed to evaluate human mobility patterns in Zambia, which to date, have largely been unexplored.Overall, mobility decreased during the first Wave (March-May 2020) following the implementation of strict NPIs including school and workplace closures, stay at home orders, and restrictions on gatherings.Research done by Gibbs et al. (2020) suggest similar conclusions of reduced mobility during 2020 with greater reductions seen at the beginning of the pandemic [19].During this time, COVID-19 cases were rising, but were maintained at lower levels compared to other countries in Africa [20].However, following the first wave, mobility returned to baseline and remained at baseline levels for the duration of 2020.Despite subsequent waves of COVID-19 and the continuation of less stringent NPIs, mobility did not drop at the same magnitudes seen during Wave 1. Further, COVID-19 cases continued to rise and fall while mobility remained at baseline levels suggesting increases in cases did not produce reductions in mobility.Rather, it is possible that NPIs played a larger role in influencing mobility, especially during Wave 1.
Estimating human mobility during the pandemic has been of large interest and companies such as Google and mobile service providers have been important sources of data.In this analysis, we highlight the differences between these two sources and conclude that mobile phone data contributes more spatially refined data by collecting data based on phone towers and not app usage, which is beneficial for mapping mobility in more rural areas.Conversely, Google Mobility Reports rely on app usage to illustrate weekly national peaks and dips categorized by type of mobility, which are not observed in the mobile phone data.Overall reductions suggested by Google Mobility Reports likely underestimate reductions in less frequently traveled routes as well as reductions seen nationally.The accessibility of these data makes it an attractive source for public health officials.However, these data may be biased; mobile phone data relies on device ownership and Google data are reliant upon device owners having their location services on and utilization of Google Maps.Additionally, Google Mobility Reports only share percent changes from baseline as opposed to point estimates for number of trips.Despite potential drawbacks, some studies have successfully utilized Google Mobility Reports to broadly evaluate the effectiveness of NPIs on human mobility and the impact on infectious disease spread [7].However, when evaluating infectious disease transmission dynamics, more spatially refined data are needed to inform tailored and effective mitigation strategies; data of which cannot be obtained from national estimates provided by Google Mobility Reports.In this analysis we demonstrate that differences in data granularity can impact the results of disease transmission modelling.Overall, both types of data have benefits and drawbacks, but both sources have been critical to aid researchers and public health officials in mapping the impact of NPIs and simulating trajectories of SARS-CoV-2 transmission.
There are a few limitations of this analysis.First, the data comes from only one of three operational mobile phone providers in Zambia and, though having nation-wide coverage, we do not know the proportion of subscribers in the country and may not be representative of the entire population.Second, we did not have data available for 2019, or pre-pandemic, time periods.In this analysis we used a period where mobility leveled out (September-December 2020) to serve as a proxy for baseline mobility, which may be an underestimation of true baseline mobility values.Lastly, due to Google Mobility Reports only providing national estimates, we were unable to extract actual values for district routes; instead, we took a proportion of the baseline data (September-December 2020) based on the national reduction estimates.
Overall, strict NPIs impacted mobility in the early stages of the pandemic.Given the risk of variants, it is important to continue the consideration of NPIs to curb COVID-19 outbreaks.Human mobility plays a large role in the propagation of infectious diseases and continues to be a key element in strategies aimed at reducing transmission.The continued efforts of companies such as Google and mobile phone services providers to contribute data are critical for public health professionals to illustrate human mobility and behavior.By mapping human mobility, more effective public health measures can be deployed to mitigate spread of infectious diseases, such as the current COVID-19 pandemic.While SARS-CoV-2 transmission is likely to transition into an endemic disease, it remains critical to consider the intersection of NPIs and human mobility as we navigate through the future of the COVID-19.

Fig 1 .
Fig 1. Overview of districts included, COVID-19 case counts, and baseline mobility.(A) The map of Zambia and bordering countries.Zambia districts and provinces (outlined in black) are labeled with districts without mobile phone data shaded in grey.(B) The time series of COVID-19 cases in Zambia from March 2020-August 2021 with the dates when key NPIs instated shown in blue bars.(C) The districts of Zambia displaying the log number of trips between districts for the top 10% most frequently traveled routes.Data were normalized to a fixed range between 0 and 1 by scaling to minimum and maximum values, and are shown for both baseline (September-December 2020) and Wave 1 (March-May 2020).https://doi.org/10.1371/journal.pgph.0000892.g001

Fig 2 .
Fig 2. Inter-province versus intra-province daily trips and return to baseline mobility levels following implementation of NPIs.(A) The average daily trips between and (B) within provinces from March 2020 to July 2021.Points represent daily average mobility estimates and solid line is the smoothed conditional average of the points with the 95% confidence interval shaded.(C) The ratio of inter-to intra-province travel calculated as the number of trips between provinces divided by the number of trips within a province.Points represent the daily ratio of inter-versus intra-province travel and the solid line represents the smoothed conditional average of the points.(D) The proportion of inter-and intra-province routes that have recovered to baseline values over time (in months).Baseline mobility and recovery time was calculated for individual routes.Overall (dashed line) represents the average across inter-and intraprovince routes.https://doi.org/10.1371/journal.pgph.0000892.g002

Fig 3 .
Fig 3. National holidays or school breaks impact on overall mobility.(A) Total number of daily trips per day (per 10k) comparing trips on days with no holiday and national holidays.(B) Total number of daily trips per day (per 10k) comparing trips on days with no holiday and school breaks.(C) Total number of daily trips per day (per 10k) comparing trips on days with no holiday and any holiday (defined as school break and/or national holiday).https://doi.org/10.1371/journal.pgph.0000892.g003 Fig 5B-5D and S8 Fig).Across all districts, Google data

Fig 5 .
Fig 5. Comparing mobile phone data and google mobility reports in overall mobility reporting and COVID-19 case predictions.(A) The percent change from baseline mobility during Wave 1 (March-May 2020) for Google Mobility reports (pink) and Mobile Phone Data (blue).Mobile Phone Data include all intra-province trips.Dashed lines represent implementation of key NPIs such as school and workplace closures.Large dips shown by Google Mobility reports are national holidays where there were large reductions in mobility to and from workplaces.The frequency of districts in which either Google or Mobile Phone data predicted a faster rise in SARS-CoV-2 cases across (B) low (�2,241 total baseline trips), (C) moderately (2,240-7,740 total baseline trips), and (D) highly trafficked routes (�7,741 total baseline trips).Each data point that makes up the range of the boxplot are the simulation runs.In other words, the y-axis "frequency" refers to the number of simulations that produced the result of any given mobility dataset predicted cases earlier.https://doi.org/10.1371/journal.pgph.0000892.g005 At baseline (averaged from 01 September-31 December 2020) a total of 880,337 average daily trips were recorded between districts (range 1-200,277 daily trips per route) by the mobile phone data.Most trips (46%) occurred within a province with trips from the capital province, Lusaka, accounting for 12% of total trips.As in other settings, travel between districts greatly decreased (40.8%,IQR: 21-70%) during the first set of national non-pharmaceutical interventions including school closures in late March (see Materials and methods, (Fig1C).98%(IQR:37-73%) relative to baseline.Similarly, the most frequently traveled district routes during baseline were also those that saw the greatest relative reductions in travel (average 70%, IQR: 68-75%) (Fig1C).
Southern: 78%, Central-Copperbelt: 90%) which corresponds to a decrease up to 2,208 daily trips per route (range 802-2,208).Less traveled routes overall experienced a lower relative reduction in travel (on average 30%, see (S7 Fig) with a few (7/51) locations also showing an increase in travel.Inter-district travel showed similar patterns to inter-province travel with reductions of 1-